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ABSTRACT 

Reflectance calculations of soybeans and com crops at two times 
during the growing season Indicate that the high sensitivity of the 
thematic mapper mld^lnfrared band to exposed bare soil between soybean 
rows is most likely responsible for early season spectral discrimination 
of corn and soybean crops by this band. 

INTRODUCTION AND SUMMARY 

Time series plots of "Greenness" and "Brightness" for corn and 
soybeans derived from the Lanasat Multispectral Scanner (MSS) bands for 
the identification of crops have previously shown that, up to about mid 
season, the time series plots are not sufficiently different to distin- 
guish soybean fields from corn fields. Roller, et. al (1980* These 
plots diverged later (after July) in the growing season. 

Corn and soybean crop reflectances for mid end late season condi- 
tions are modeled using a row crop reflectance model, Suits (1983), in 
order to determine the cause for this similarity in early season band 
spectra and to assess the potential of Thematic Mapper (TM) bands for 
distinguishing these two crops at mid season. 



The reflectance model analysis for mid July shows that the MSS mean 
vectors for soybeans and corn are nearly col linear. A small variation in 
the soybean foliage extension into the inter-row region will result in 
Greenness and Brightness values quite close to the mean values for corn. 
The mean vectors for corn and soybeans change direction after July. The 
Greenness of corn remains about the same but soybean Greenness increases 
significantly. The analysis suggests that the divergence is due to the 
difference in the maturing and scenesing features of the two crops. 

A two-band plot using the 0.79 to 0.90 and the 1.55 to 1.75 ym bands 
of the Thematic Mapper shows that the mid July TM mean vectors are not 
coll inear but are separated by 16 degrees, offering the potential for 
early discrimination. The band reflectance for the MSS and TM near- 
infrared bands are nearly the same. The difference is in the sensitivity 
of the 1.55 to 1.75 ym band of the TM to exposed soil in soybeans. 

The analysis begins w^th the collection of crop descriptors for corn 
and soybean fields for mid July and mid August. These descriptors are 
used as input to the row crop reflectance model and the reflectances are 
calculated under sun angles appropriate to a midwest-U.S. satellite 
over-pass time. Two-band plots show the location in spectral space of 
the crop mean vectors for early soybean, early corn, late soybean and 
late corn. 

A re'Hectance sensitivity calculation yields the expected shift in 
the soybean mean vector for hypothetical field feature variations. The 
feature variation which causes the largest shift towards or away from the 
corn mean vector is used to identify the feature which is primarily 
responsible for the spectral separation of the two crops. 


FIELD DESCRIPTORS 

Two categories of field descriptors are required fo*' the row crop 
reflectance model Input — spectral properties and geometric properties. 

The hemispherical spectral reflectance and transmittance of the 
plant components — green leaves, cnlorotlc leaves, necrotic leaves, and 
stems — were obtained by laboratory measurements on samples supplied by 
G. Saflr at Michigan State University. Soil reflectance was obtained 
similarly. The soil reflectance was a Type 1 spectrum by Condi t (1970) 
classification. 

The geometric data — the mean component projected areas (horizontal 
and vertical), and the number of such components per unit volume of crop 
canopy (number density) — was obtained from a combination of field 
observations and reported information. 

Leaf orientation for corn was obtained from field observations at 
Michigan State University. Soybean leaf orientation was derived by G. 
Safir from Blad (1970). The canopy height, row spacing, green leaf area 
index (LAI), and canopy layer thicknesses were derived from interpreta- 
tion of LARS (1980) field photographs and field reports for the two times 
of interest during the growing season, mid July and mid August. 

Along with these descriptors, an addititonal description of the row 
structure of the crop is needed. The number of components per unit 
canopy volume is specified to be the field average. Across rows, the 
number density of components changes from maximum on the row to some 
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reduced value between rows. A simple linear variation of number density 
was chosen for corn where the number density of components are reduced 
linearly with distance towards the adjacent row. The variable, foliage 
extension (FE), is the fraction of row spacing at which the nuirtfjer 
density of components extending from one row to the next becomes zero. 
Each row is assumed to behave symmetrically so that at FE * 0.5, the 
foliage from one row progressively becomes more sparse and is zero at the 
point where the foliage from the opposite row begins to increase from 
zero density towards its row center. A foliage extension of 1.0 requires 
that the number density of one row becomes zero at the opposite row 
center. The sum of the number densities between rows due to extension 
fron adjacent rows is constant and the field is uniform. 

The soybean foliage is distributed differently. The LARS field 
photographs indicate that the distribution of components across rows is 
more nearly in the fashion of a rectangular distribution where foliage 
extending into the inter-row regions appears to be uniform up to the 
point where number density falls to zero abruptly. For a soybean field 
at FE = 0.5, the field is uniform with foliage from opposite rows meeting 
in the middle. At FE < 0.5, a clear strip of bare soil can be seen 
between soybean rows. 

The basic parameters for the early corn and soybeans which were used 
in the analysis are shown in the Appendix. Changes only in geometric 
parameters were made in order to model late corn and soybeans. The LARS 
field photographs indicated that the late corn canopy differs from the 
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*iar1y corn canopy primarily In the lower canopy layer where after mid 
July leaf senescence begins progressively upward towards the top of the 
canopy. The top layer green leaf area Index appears to be constant. 
Consequently, the July layer thicknesses of 0,25 m (lower) and 1.75 m 
(upper) were changed to 1,0 m (lower) and 1.0 m (upper). The lower layer 
contains the mixture of green, chlorotic, and necrotic leaves in the leaf 
slough-off region. 

The LARS data shows that the late soybean canopy shows no evidence 
of senescence during that time but instead soybeans have increased their 
green leaf area index from 3.5 to 5 and have Increased the foliage 
extension to FE » 0,5 to form a unifonn canopy. Thus, the late soybean 
field parameters are: green leaf projected areas are Increased by 5/3.5 

and FE « 0.5. 

ANALYSIS OF BAND SPECTRA 

The calculation of the reflectance spectra, p(x) was accomplished 
using the row crop reflectance model. Suits (1983), from 0.4 to 0.8 ym in 
steps of 0.01 ym, and from 0.8 to 1.1 ym, from 1.5 to 1.75 ym, from 2.0 
to 2.35 in steps of 0.05 ym. Band reflectances were calculated from 
these spectra using the relation, 

^2 

p(band) » ( 1 /(X 2 - X^)) J' p{\) dX 

^1 

where x^ and X 2 are the nominal band limits. 
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The four MSS bands and the six daylight reflectance bands of TM were 
used. For convenience, these bands are renund>ered and assigned nominal 
band limits as shown In Table 1. 

The resulting band spectra for the four cases — early corn (Cl), 
early soybeans (SI), late corn (C2), and late soybeans (S2) — are shown 
In Table 2. If the field parameters describing the fields wore the 
iverage values for a hypothetical population of similar fields, then 
these band spectra could be the mean vectors for the respective popula- 
tions. A two-band spectral subspace was selected In order to show the 
relationships graphically between these hypothetical crop mean vectors. 
The subspace. Band 2 and Band 4, was chosen for the MSS and the subspace. 
Band 8 and Band 9, was chosen for the TM. Figure 1 shows the MSS two- 
band mean vectors for the four fields. The scale for Band 2, the red 

band, was greatly expanded for clarity. Notice that the mean vectors for 

Cl and SI are nearly col 11 near. Using the dot product between Cl and SI 
vectors, one can determine the angle between to be about 0.3 degrees. 

Each mean vector points to the center of its respective hypothetical 
scatter plot representing the results for the respective hypothetical 
populations. 

Although the actual shape of the cluster of points which would be 
made by a scatter plot Is not obtained from the calculation, experience 

with both Landsat and airborne scanner scatter plot data has shown that 

the best fitting ellipse around the clusters (assuming a Gaussian 
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TABLE 1. BAND NUMBERS AND BAND LIMITS FOR THIS ANALYSIS 


Band 

Nun^er 

Band 

Name 

Allas 

ym 

yin 

1 

MSS 4 

green band 

0.50 

0.60 

2 

MSS 5 

red band 

0.60 

0.70 

3 

MSS 6 

IR 1 

0.70 

0.80 

4 

MSS 7 

IR 2 

0.80 

1.00 

5 

TM 1 

bl ue-green 

0.45 

0.52 

6 

TM 2 

green band 

0.53 

0.61 

7 

TM 3 

red band 

0.63 

0.69 

8 

TM 4 

near IR 

0.79 

0.90 

9 

TM 5 

m1d IR 

1.55 

1.75 

10 

TM 7 


2.10 

2.35 
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TABLE 2. CORN AND SOYBEAN MEAN VECTORS (X) 
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BAND 4 (X) 



BAND 2 (%) 

FIGURE 1. CORN AND SOYBEAN MEAN VECTORS IN THE MSS {2, 4) SPECTRAL SUBSPACE. 
Cl * early corn; C2 * late corn; SI * early soybeans, 

S2 ■ late soybeans. The angle between SI and Cl Is 0.3® 
with Euclidian reflectance distance of 6.2^. The angle 
between S2 and C2 is 3* with Euclidian reflectance 
distance of 12.5%. (Notice the expanded scale for 
Band 2.) 


distributions) Is almost always oriented In such a way that the major 
axis of the ellipse Is near the direction of the mean vector. The over- 
lap of the scatter plots win result In errors In classification when 
multlspectral discrimination algorithms are applied. The case of Cl and 
SI would probably produce greatly overlapping scatter plots In the major 
axis direction. 

The late season case, C2 and S2, shows that the mean vectors have 
separated both In distance and In direction. Although the angle between 
has Increased to only 3 degrees, yet, because of the probable shape of 
the scatter plots, the change In angle may be the most Influential In 
making spectral discrimination algorithms more accurate for late corn and 
soybeans. 

Figure 2 shows the same four fields using the TM two-band subspace. 
The band spectra scales are the same in this case. Notice the wide 
separation for the early fields. Cl and SI. The angle between mean 
vectors is about 16 degrees. 

Such a large difference in direction cannot, by Itself, lead to the 
conclusion that early season discrimination would be Improved because the 
minor axes of the scatter plots may also be larger. However, It Is 
certainly a change in the right direction. The results of experiments by 
Unger and Goward (1982) showed that, using Band 9 Instead of the MSS 
bands, early spectral separation of corn and soybeans was possible In 
their sample. The question remains as to the physical cause for the 
early separation. 
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BAND 8 {%) 



FIGURE 2. CORN AND SOYBEAN MEAN VECTORS IN THE TM (8. 9) SPECTRAL SUBSPACE. 
Cl * early corn, C2 ■ late corn; SI * early soybeans, 

S2 » late soybeans. The angle between SI and Cl is 16.1® 
with Euclidian reflectance distance of 12.4%. The angle 
between S2 and C2 is 5.4° 'with Euclidian reflectance 
distance of 17.6%. 
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Could the cause be accidental or coincidental to the partlcut.r circum- 
stances of the Unger and Qoward field measurements or to the particular 
choices of field parameters In this calculation? 

SENSITIVITY ANALYSIS 

If one could identify the physical feature of these fields which Is 
primarily responsible for the separation using the TM sub-bands, then one 
could assess the possibility that the cause may or may not be purely 
circumstantial for particular choices of test area or choices of modeling 
parameters. A sensitivity analysis was done In order to Identify this 
feature. 

The mean vectors were recalculated with certain field features 
changed one at a time. The nevy mean vectors, p” , are related to the 
previous mean vectors, p, by 


*p' * "p + S^J) 

where I(J) ■ shift vector caused by a change in feature J. 

Changing some field features required only a change in one of the 
model parameters; others required a combined change. 
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Ten field features were chosen for variation: 

1. Percent change in plant density (%PD). 

2. Percent change in green LAI (%LAI). 

3. Degrees change in leaf slope (LS). 

4. Percent change in soil reflectance (%RH0). 

5. Change in foliage extension (aFE). 

m 

6. Percent change to chlorosis (XCHL). 

7. Percent change to necrosis (%NEC). 

8. Percent change in row spacing 

9. Degrees change in row azimuth (aAZ). 

10. Percent change in canopy height (XHCAN). 

With the exception of changes in chlorotic and necrotic leaf areas, 
individual features were made to change with all other features held 
fixed. For example, a change in row spacing will not alter the foliage 
extension. The foliage would be more thinly distributed than before but 
would reach the same fractional distance across rows. The percent change 
to chlorosis ana necrosis requires that green leaf area be turned into 
chlorotic or necrotic leaf area. Green LAI is reduced and is transferred 
as increments to the chlorotic or necrotic leaf area without any other 
change. 

The changes in field features which were used and the shift vectors 
for each change are shown In Table 3 for early soybeans. 
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One may observe that Band 9 In Table 3 appears to respond to a soil 
reflectance Increment, Feature 4, more than to any other field feature 
change. Also Band 9 responds to soil reflectance with greater sensiti- 
vity to that feature than does any other band. The response of Band 10 
is a close second. Further, one may observe that a change in row azinwth 
from east-west to north-south causes a response in the opposite direction 
in both Bands 9 and 10 with Band 9 being second to Band 10. The sun azi- 
muth of 104* and poUr angle of 41® causes some intcr-row soil illumin- 
ation. A change to north-south rows places the sun at approximately the 
cross-row direction and should shadow much of the inter-row soil. The 
shift due to a change of row direction from east-west to north-south will 
tend to shift the early soybeans mean vector in the direction of the 
early corn mean vector. 

One may conclude that the primary cause for the separation of early 
corn and soybean mean vectors using the TM subspace bands is the high 
sensitivity of Band 9 to tne exposed and illuminated soil in the soybean 
field. The MSS subspace Bands 2 and 4 also respond to the same feature 
and in the same sense but the red band is less sensitive than the TM 
Band 9. 

The sensitivity of the two TM lor.n wavelength bands to soil varia- 
tions has been previously noted by Crist (1983). He reports that, where 
Kauth and Thomas (1976) described a "plane of soils" in the MSS Tassel ed 
Cap space, in MSS data virtually all the variation normally encountered 
in soils was confined to the Brightness direction so as to define a "soil 
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TABLE 3. EARLY SOYBEAN SHIFT VECTORS {%) 
BAND 



15 


line" Instead of a plane. In the analogous TM Tassel ed Cap space, a true 
"plane of soils" can be seen which Is due to the contribution of the two 
TM long wavelength bands. 

This sensitivity to exposed soils Is reduced when Illumination fails 
to reach the soil so that the separation of early soybeans and corn can 
be reduced in cases where solar elevation Is both low and the sun Is near 
the across-row direction at overpass time. The separation will then be 
less pronounced. 

acknowledgements 

This work was sponsored by NASA/JSC Earth Resources Research 
Division. The author wishes to thank Prof. G. Safir and Dr. W. Mai 11 a 
for their assistance In supplying data and data sources. 


16 



APPENDIX 


FIELD PARAMETERS 


The parameters used for the early mean fields are shown in Table 4. 
The sun is at polar angle, 41®, and azimuth from north, 104®, to repre- 
sent a 10 AM EST satellite overpass in the midwest in mid July. The view 
direction is the nadir. 
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TABLE 4. FIELD PARAMETERS FOR EARLY CORN 
AND SOYBEANS (MID JULY) 


Parameter 

CORN 

SOYBEANS 

Planting Date 

Mid May 

Mid May 

Plant Density 

5.0 m’^ 

27.4 m"^ 

Row Direction 

East-West 

East-West 

Row spacing 

0.75 m 

0.75 m 

FE 

0.6 (triangular) 

0.3125 (rectangular) 

Canopy Height 

2.0 m 

0.45 m 

Top Layer 

1.75 m 

0.30 m 

Bottom Layer 

0.25 m 

0.15 m 

Green LAI 

3.5 

3.5 

Layer 1 Slope 

60® 

52® 

Layer 2 Slope 

75® 

47.4® 

Extinction 

Coefficients 

Top Layer 

Green Leaf H 

0.874 m"-^ 

4.80 m"^ 

V 

1.514 m“^ 

6.13 m"^ 

Chlorotic H 

0.0875 m‘^ 

0.001 m*^ 

V 

0.154 m’^ 

0.001 m'^ 

Necrotic H 

0.0875 m’^ 

0.001 m"^ 

V 

0.154 m’^ 

0.001 m‘^ 

Stem H 

0.0 

0.049 m"^ 

V 

0.303 m“^ 

0.135 m‘^ 

Bottom Layer 

Green Leaf H 

0.456 m'-^ 

5.27 m"^ 

V 

1.700 m‘^ 

5.73 m’^ 

Chlorotic H 

0.080 m“^ 

0.02 m'^ 

V 

0.348 m"^ 

0.02 m'^ 

Necrotic H 

0.080 m'^ 

0.02 m‘^ 

V 

0.748 m"^ 

0.02 m'^ 

Stem H 

0.0 

0.025 m“^ 

V 

0.280 m‘^ 

0.068 m’^ 

Soil 

Condi t Type 1 

Condit Type 1 
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